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In this work we investigate the dynamics of random walk processes on scale- free networks in a short 
to moderate time scale. We perform extensive simulations for the calculation of the mean squared 
displacement, the network coverage and the survival probability on a network with a concentration 
c of static traps. We show that the random walkers remain close to their origin, but cover a large 
part of the network at the same time. This behavior is markedly different than usual random 
walk processes in the literature. For the trapping problem we numerically compute <E>(n, c), the 
survival probability of mobile species at time n, as a function of the concentration of trap nodes, 
c. Comparison of our results to the Rosenstock approximation indicate that this is an adequate 
description for networks with 2 < 7 < 3 and yield an exponential decay. For 7 > 3 the behavior is 
more complicated and one needs to employ a truncated cumulant expansion. 



PACS numbers: 05.40.Fb, 82.20.-w, 89. 75. Da 



I. INTRODUCTION 

Scale-free networks have been widely studied during 
the recent years, mainly because of their connection to 
a plethora of real- world structures 0, 01 ■ These net- 
works are made by nodes connected to each other via 
links, which may be directed or undirected (in the present 
work we only deal with undirected Hnks). Studies of their 
structure show that most of them possess the small- world 
property, i.e. the mean path length is extremely small 
and every node can be reached by following a path con- 
sisting from a very small number of nodes, as compared 
to the case of lattice systems. A special feature, though, 
that distinguishes this class of networks is the fact that 
the probability distribution for a node to have k links to 
other nodes obeys a power law: 



P{k) 



(1) 



where 7 is a parameter that measures how densely con- 
nected a network can be. There is a wide range of real- 
life networks yj, |2] that have been shown to follow this 
power-law form in their connectivity, including networks 
in nature, such as the cell, metabolic networks and the 
food web, artificial networks such as the Internet, the 
WWW and power grids, or even social networks, such as 
sexual partnership networks. 

The scale-free networks, termed after the absence of 
a characteristic typical node connectivity, exhibit many 
unusual properties as compared to simple lattice models, 
random graphs, or even small- world (Watts-Strogatz Q) 
networks. This scale-free character results in the exis- 
tence of a small number of nodes which are connected to 
a large number of other nodes. These super-connected 
nodes (termed 'hubs') have been shown to have a central 
place in the interpretation of many of the network prop- 
erties. A lot of work has been devoted in the literature to 
the study of static properties of the networks, while inter- 
est is growing for dynamical properties on these networks. 
Recently U , we presented results for the absence of ki- 
netic effects in reaction-diffusion processes taking place 



on scale-free networks. In this work, we study a num- 
ber of random walk properties, including mean-squared 
displacement, network coverage, and trapping processes 
on scale-free networks of varying connectivities. Trap- 
ping has been considered in the past as a model for en- 
ergy transfer, but also in a more general frame in the 
context of networks, as a model for the probability of 
reaching targets located on the network in a given con- 
centration via random moves iS'. It is of interest, thus, 
to study the mechanism, the effects of connectivity, con- 
centration, size, etc on such structures that exhibit these 
unique properties. Our results refer to small to moderate 
time regimes, where we are still far from the asymptotic 
limit. This limit has been known to be very hard to reach 
in regular lattices, too, and cannot be predicted by direct 
simulation techniques. 



II. RANDOM WALKS AND THE TRAPPING 
PROBLEM 



One of the most basic quantities in the random walk 
theory is the mean-squared displacement {R^{n)) of a 
particle diffusing in a given space, which is a measure 
of the distance R covered by a typical random walker 
after performing n steps. In most cases, this quantity is 
described by an expression of the form: 



(i?2(n)) - n'' 



(2) 



The value of the parameter a classifies the type of dif- 
fusion into normal linear diffusion (a = 1), sub-diffusion 
(a < 1), or super-linear diffusion (a > 1). Of course, 
when we consider distinct time steps and nearest neigh- 
bor lattice hops the maximum value of a can be 2, i.e. 
a completely biased walk where the random walker con- 
tinuously moves away from the origin. Recently 6], the 
mean squared displacement was studied in small world 
networks, where it was shown that diffusion is linear and 
results were found to collapse under a proposed scaling. 



The behavior of a random walk is also characterized 
by the coverage of the space, as expressed by the average 
number of distinct sites visited (Sn) after n steps. In 
regular Euclidean lattices this quantity follows a power 
law with the number of steps, except in the case of two 
dimensions, where logarithmic corrections appear in the 
denominator {{Sn) ^ n/lnn). In one dimension (5„) ^ 
y/n, and in dimensions higher than two (S'„) ^ n, and the 
number of sites visited grows linearly with the number 
of steps n, since the random walker practically visits at 
each step a new site. In infinite dimensions, of course, 
the number (5„) of visited sites is equal to n, since there 
is no revisitation of sites during the walk, and the walker 
covers the largest possible area. In small-world networks 
a scaling ansatz was proposed 6], which was verified by 
simulations, and {Sn) shows a transition from a slope 0.5 
(one-dimensional behavior) to a slope 1 (d > 2 behavior). 

An important process related to random walk theory 
is trapping |3, l3l ■ Trapping reactions have been widely 
studied in the frame of physical chemistry, as part of 
the general reaction-diffusion scheme. The general idea 
includes two different species A and B, which diffuse 
freely in a given space and upon proximity they react ac- 
cording to A-l-B^-B. Many different variations describe 
a plethora of physical phenomena. In this paper we deal 
with the special case of the trapping problem where B 
particles are immobile. The simplest mean-field ana- 
lytical treatment predicts a simple exponential decay in 
the density of A's, while the earlier contributions to the 
subject go back to Smoluchowski 9], who was the first 
to attempt to relate the macroscopic behavior with the 
microscopic picture by taking into account local densit 
fluctuations. However, over the years a lot of work |3, 
has been devoted to the trapping problem which, even in 
its simplest form, was shown to yield a rich diversity of 
results, with varying behavior over different geometries, 
dimensionalities and time regimes. 

The main property monitored during such a process is 
the survival probability <I>(n, c), which denotes the prob- 
ability that a particle A survives after performing n steps 
in a space which includes traps B with a concentration 
c. It is well-known that $ behaves differently in dif- 
ferent dimensions, as well as in different time-regimes. 
The problem was studied in regular lattices and in frac- 
tal spaces [3, Q, and, recently, in small- world networks by 
Blumen and Jasch [l3,|ll[i3. 

The simplest treatment of the trapping problem on 
a lattice assumes that when a random walker has per- 
formed n steps and has visited Sn different lattice sites 
at least once, the probability that it has not yet been 
trapped is equal to (1 — c)^" , where c is the trap concen- 
tration. When this quantity is averaged over all different 
possible walks and trap configurations the resulting sur- 
vival probability will be equal to 



tutes the above expression with the typical value of the 
distribution, i.e. 



$(n,c) = ((l-c)^") = (e-^^") 



(3) 



where A = — ln(l — c). A simplification of this equation 
was first proposed by Rosenstock [ij] and simply substi- 



$(n,c) 



-A(S„) 



(4) 



This approximation has the advantage that the mean 
value of the number of sites visited {Sn) is well known [lj| 
for practically all dimensionalities (including e.g. fractal 
ones). Notice that the Rosenstock approximation does 
not necessarily imply a simple exponential decay, except 
in the case where {Sn) ^ n. The formula predicts sim- 
ple exponential decay of the survival probability with the 
number of steps n only for d > 3, and exponential de- 
pendence on ^/n in d = 1. In 2 dimensions the predicted 
behavior is rather complex, with logarithmic corrections 
in the exponent. The applicability of Eq. Q is limited 
to short-times and/or not too large trap concentrations. 
When the survival probability becomes low enough, this 
expression deviates significantly from the correct behav- 
ior. 

A significant improvement was possible by the use of 
averaged quantities, known as cumulants, where the av- 
eraged quantity of Eq. (Q can be written as a function 
of the cumulant generating function [13 : 



Kj{\n)^Y.^- 



Ayjk,{n) 



(5) 



where kj{n) are the cumulants, which are associated to 
the moments of S'„, e.g. fci(n) — {Sn), k2{n) = {Sn) — 
{Sn)^, etc. The expression Q for the survival probability 
then simply becomes 



^jin,c) = exp{Kj{X,n)) 



(6) 



Improved accuracy can be obtained upon increasing the 
truncation order J. In theory, the knowledge of all the 
moments (J -^ 00) for the Sn distribution is required for 
the use of JS)). These moments are known analytically 
only in one-dimensional lattices, while for d > 1 usually 
the first 2-4 moments are used. 

A detailed analytical treatment of the problem was 
performed by Donsker and Varadhan •Id'l , who were able 
to predict the asymptotic behavior of the survival prob- 
ability as 



lim ^{n,c) — ex.p{~KdX^+''n''+^) . 



(7) 



The positive constant Kd depends only on the dimension- 
ality and the structure of the lattice. This asymptotic ex- 
pression does not provide any information on when the 
asymptotic limit is reached. Since it has been observed 
that the Rosenstock approximation describes quite well 
the high-$ regime, it is obvious that a crossover to the 
Donsker- Varadhan behavior will take place. The location 
of this crossover has been studied in detail [13, Hsl u^ , 
and it was shown that only with indirect methods it is 
possible to reach the asymptotic limit. 



This asymptotic behavior has also been explained via 
heuristic arguments. The slow relaxation of $ at long 
times is due to an interplay of two different factors. First, 
mean-field treatments assume a uniform trap distribution 
over the entire space. This is not strictly true, though, 
and for large enough sizes it is possible to find very ex- 
tended trap-free regions. A random walker in such a 
region will survive for extremely long times compared to 
walkers in normal regions and will thus determine the 
asymptotic behavior. The second factor is due to unusu- 
ally 'compact' random walks, which revisit many times 
the same sites, and thus result to a very small value of 
S„, even at longer times. 

Recently, a number of papers were published concern- 
ing trapping on a version of the small-world networks 
[Hi mi Il2 | . These networks, first proposed as a model by 
Watts and Strogatz J3], are one-dimensional rings where 
additional links are inserted between two random sites 
with a given probability. It was shown that the results 
represent a fine interplay between pure order and pure 
disorder statistics. Initially, the walkers feel only the 
presence of the one-dimensional lattice, but at longer 
times the behavior of the survival probability follows that 
of an open tree structure. The decays of the survival 
probability were clearly not exponential, and the cumu- 
lants description did not yield accurate coincidence with 
the numerical results in all of the studied cases. 

In this work, we extend the above mentioned random 
walk problems (mean-squared displacement, coverage, 
and trapping) in the case where the underlying struc- 
ture is a scale-free network, obeying a power-law in the 
nodes connectivity distribution. The random walkers are 
located on the nodes and can only move along the links of 
this network. In the case of the trapping model certain 
nodes are designated as traps, having a concentration 
c. We present computer simulations results for different 
network connectivities and compare them to the known 
lattice behavior. 



connectivity. However, it is possible that at the last 
stages of the construction we will reach a dead-end where 
no further links may be established according to the 
above rules. In this case we simply ignore the 'hands' 
that cannot be connected, since their number is always 
very small and the structure of the network is not influ- 
enced at all. 

The largest cluster in the network is identified via the 
use of a spreading algorithm. We start with a random 
node and mark it with a label, say X. We then mark all 
the nodes connected to this node by X, and proceed it- 
eratively by labeling their neighbors, etc, until the whole 
cluster has been labeled. We then choose another ran- 
dom node that has not been labeled by X, which means 
that it belongs to a different cluster. We mark it by Y 
and again spread this labeling throughout this cluster. 
When the entire network has been labeled we can easily 
identify the largest cluster from its size. 

All random walks in this paper take place on the largest 
cluster of the network via the following algorithm. We 
place a random walker on a randomly chosen site i of 
the largest cluster. This site has a connectivity fc^. At 
each Monte-Carlo step the walker makes a jump towards 
a node connected to i (i.e. nearest neighbor) with prob- 
ability 1/ki. This process gives a Markovian walk, since 
each step is independent of all previous steps. Distances 
on the network are measured according to the shortest- 
length path between two nodes, and the displacement R 
of a walker is calculated relative to the initial point. 

For the trapping problem, we randomly choose a per- 
centage c of the network nodes and designate them as 
traps. A random walker is placed on a random non-trap 
node and performs the procedure described above until 
it meets a trap. In this case, it is annihilated and the 
time n to trapping is recorded. We repeat the same pro- 
cess for many independent random walkers and different 
networks, and we construct a histogram iJ(n, c) of the 
number of walks that last exactly n steps. Then, the 
survival probability is simply given by 



III. THE MODEL 

The construction of a scale-free network follows the 
MoUoy-Reed scheme |23- First, we fix the number of 
nodes N in the system and the 7 parameter, character- 
istic of the particular network. By using the transforma- 
tion method we select N random numbers from the k^'^ 
distribution, so that each node i is assigned a number of 
links ki from the above distribution. The value of k lies 
in the range from /cmin = 1 (lower cutoff) to fcmax = iV— 1 
(no upper cutoff value is used for k). 

Initially, no links are established in the system. Each 
node i extends ki 'hands' towards all other nodes. We 
randomly select two such 'hands' (that do not belong in 
the same node) and connect them creating thus a link. 
No double links are allowed, so that if two nodes are al- 
ready connected this link is rejected. We continue this 
process until all nodes have reached their pre-assigned 



1 
$(n, c) ^ 1 - — ^ i/(TO, c) , 

m— 1 



(8) 



where AI is the total number of independent random 
walks sampled. 

Typically, 100 different networks with N = 10^ nodes 
were created, and 10^-10* different random origins were 
chosen on each network. Thus, results were averaged over 
10^-10^ diiferent realizations of the walk. 



IV. RESULTS 

A. Mean squared displacement 

We first study the mean squared displacement (i?^) 
of a random walker on a scale-free network. Since the 
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FIG. 1: Mean squared displacement as a function of time for 
networks of 7 = 2.0, 2.5, 3.0 and 3.5 (shown on figure). The 
network size in all cases is A'' = 10^ nodes, and for 7 = 3.0 we 
also present results for networks of size A*' — 10'* and A'^ — 10^ 
(bottom to top). The dotted lines represent slopes of 1 and 
2. Inset: Normalized {R'^{n)) curves, so that asymptotically 
all curves converge to 1, for different 7 values (shown on the 
plot). Dotted lines represent best-fit lines with slopes (left to 
right) 1.8, 1.5, 1.1, and 0.9. 



networks that we study are not embedded in a regu- 
lar Euclidean space, this quantity does not measure how 
far in Euclidean space the walker travels, but rather the 
minimum number of hops needed in order to return to 
its origin. The first important feature of Fig. ^ where 
we present (i?^) as a function of time for networks of 
varying connectivity distributions, is the fact that (i?^) 
equilibrates after a few steps to a constant displacement 
value. This is a simple manifestation of the very small 
diameter of these networks, which has been shown to be 
of the order ln(ln TV) [21| . In practice, this means that 
one node can be reached from all other nodes in the net- 
work within only a few steps and the maximum possible 
distance in the network is very small compared to the 
network size. Note also that the plateau value increases 
as we increase 7, since a network which is less connected 
exhibits a larger diameter. The existence of the plateau 
is, of course, a finite-size effect. However, the size depen- 
dence is not strong, as can be seen in the figure, where 
we present results for networks with 7 = 3.0 and size 
iV = 10"*, 10^, and 10^. Although we increase the size 
of the networks by two orders of magnitude, the value 
of {R'^) increases from roughly 70 to 110, i.e. the dis- 
tance R increases almost linearly from 8 to 10. This 
logarithmic dependence shows that for all practical ap- 
plications the plateau will be present. For example, it 
has been observed [22| that the diameter of WWW (of 
size iV = 8 X 10^ and 7 = 2.45) is only 18, so that even on 
such large networks maximum distances remain small. 

In the figure inset we have rescaled the (i?^(n)) data so 
that all curves are normalized to an asymptotic value of 1 . 



FIG. 2: Number of distinct sites visited {Sn) after n steps on 
scale-free networks with (solid lines, top to bottom) 7=2.0, 
2.5, 3.0 and 3.5. The dashed line is the infinite-dimensional 
case of (Sn) — n + 1. The network size was A'^ — 10^. 



It is shown that upon varying the value of 7, diffusion on 
scale- free networks may range from superlinear to sublin- 
ear diffusion. For networks of low 7, diffusion is greatly 
enhanced. Thus, for 7 = 2.0 the walkers move away from 
the origin rapidly and the slope of (i?^) reaches a value 
of about 1.8. After only a few steps, though, the value of 
(i?^) saturates, due to the phenomenon described above. 
As we increase the value of 7 the slope of the curves 
decreases. Diffusion at early steps remains super-linear, 
until we reach a value of 7 around 3.0 where the slope 
becomes roughly equal to 1. This linear diffusion turns 
slowly into sub-linear as we further increase 7 and for 
7 = 3.5 the slope is equal to 0.9. 



B. Network coverage 

The coverage of the network by random walks is found 
to be a very efficient process. Numerical results of (Sn) 
on scale-free networks are presented in Fig. |21 The num- 
ber of sites visited increases initially with a slow rate, but 
after a crossover value the increase is almost linear. This 
asymptotic linearity is observed in all 7 values, while the 
crossover point shifts towards longer times with increas- 
ing 7. The early time slope means that the walkers ini- 
tially spend some time exploring the neighborhood they 
were created in and visit the same sites. After the first 
few steps (the exact number depends on the connectiv- 
ity of the network) they escape their initial territory and 
diffuse around the entire network. Thus, it is possible to 
continuously visit new sites, which results in the linear 
increase of (Sn)- In reference [g, (Sn) on small world 
networks was found to scale at early times with ^/n and 
asymptotically with n. In the case of scale-free networks, 
the early-time behavior is not consistent with a ^/n law, 
which would be an indication of one-dimensional behav- 



ior. For each 7 value the local environment is different, 
and this is exhibited in the different evolution of (5„) 
for low n values. The crossover, also, is located at much 
earlier times (of the order of tens of steps) as compared 
to thousands of steps which is the case reported in [g] for 
small world networks. 

The size of the network used (N =10^) was two orders 
of magnitude larger than the number of steps performed, 
in order to avoid finite size effects. Despite of this precau- 
tion, the curve of 7 = 2 seems to deviate from linearity at 
longer times. This phenomenon means that revisitation 
already starts to exhibit itself for the finite network we 
study. 

The linear growth of (Sn) is similar to the behavior 
exhibited in dendrimer structures, modeled by Cay ley 
trees. These are open structures, with every node having 
a fixed number k of connected nodes which are always 
directed away from the central core. It was also shown in 
that case [23 that (Sn) had a linear increase after a short 
early-time sublinear regime, due to the same reasons as 
here. 

Consider, now, a regular lattice that can be embedded 
in a finite d-dimensional space. In this case, it is well- 
known that unconstrained diffusion causes the random 
walker to spread in the available space, increasing both 
(i?^) and {Sn) with time. Thus, the walkers tend to in- 
crease their distance from the origin and cover new terri- 
tory. If for some reason we restrict diffusion of the walkers 
within a finite distance Re from the origin, so that (i?^) 
saturates, then the area covered will soon saturate, too, 
to a value of order R^. Diffusion on scale-free networks, 
though, is different. Although the walker is always close 
to the origin and restricted within a distance equal to the 
network diameter, new territory is continuously sampled. 
This peculiar behavior can be attributed to the existence 
of the hubs. If we consider an extreme case of a hub, 
that of a star node where all the nodes are connected 
only to the hub, then the displacement will be at most 
2 steps away, but due to the large number of nodes in 
the system the walker will be redirected to a non- visited 
site with a revisitation probability n/N, which for large 
enough systems and early to moderate times is close to 
0. The particular case of Levy flights |2j (which usu- 
ally results in enhanced diffusion with a > 1 in Eq. ^ 
can be considered similar to the process we study in this 
paper. In a Levy flight the length of a jump follows a 
power law dependence. In practice, a walker samples an 
area for a certain amount of time before performing a 
long-range jump. This jump allows then the walker to 
sample a new space. Although the areas of space visited 
can be regarded as the hubs of the present problem, the 
main difference is the displacement of the walk. In the 
case of Levy flights the 'hubs' of the system are distant 
in space among each other, while in our case the hubs are 
very close to each other, and can even be directly con- 
nected. As we have already seen, thus, the mean squared 
displacement on scale-free networks (even of large size) 
is restricted to small distances, whereas (i?^) increases 
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FIG. 3: Survival probability as a function of time, for a 
network of 7 = 2.5 and trap concentration c — 0.01. From 
right to left, the number of nodes in the network is A'^ = 10"^, 



monotonically in the case of Levy flights. 



C. Trapping 

For the trapping problem, we first examine the depen- 
dence of the survival probability $ on the system size N. 
As it can be seen in Fig. 13 for 7 = 2.5, larger networks 
yield a significantly lower survival probability. This is 
due to the higher probability of finding a node with very 
high connectivity, which is linked directly with the largest 
part of the network. Due to the power-law dependence 
the appearance of these nodes increases as we increase 
the network size. However, we can see that the $-curves 
for the larger networks {N > 10^) practically coincide. 
Moreover, this TV-dependence is much weaker for net- 
works with higher 7 values. As we can see in the plot, 
the large-size network behavior in this case is very close 
to a simple exponential decay, while smaller networks de- 
viate from this behavior. 

In Fig. ^a) we present the survival probability on the 
largest cluster for different trap concentrations as a func- 
tion of time, for networks with 7 = 2.5. For a relatively 
high trap concentration, e.g. c = 0.05, we can see that 
$ falls very rapidly and during the first 100 steps only 
a small percentage of the walkers has survived. The de- 
cay retains for the largest part an exponential character. 
In order to test the validity of the Rosenstock approx- 
imation for scale-free networks, we used the numerical 
data for (Sn) presented in Fig. [51 and computed the sur- 
vival probability <I> using Eq. Q). The results in Fig. ^ 
show that there is almost complete coincidence between 
this approximation and the simulation data. As we have 
mentioned above, the Rosenstock approximation is valid 
when mean-field features are present, and fluctuations in 
the area covered are not important. Thus, a high trap 



*(n,c) 



*(n,c) 




2000 4000 6000 8000 10000 

Time (MC steps) 



FIG. 4: Survival probability of random walkers as a func- 
tion of time, for a network of size A'^ = 10 and (a) 7 — 2.5, 
(b) 7 = 3.5. Symbols represent direct trapping simulations. 
Solid lines represent the Rosenstock approximation based on 
the (Sn) data of figure |^ Dashed lines are the results of the 
cumulant approximation, with the truncation order J indi- 
cated on the plot. From left to right, the trap concentrations 
are c = 0.05, 0.01, 0.005, and 0.001. 



concentration implies that there will be no large trap- 
free regions, since a walker can easily escape any part of 
the system. However, in the case of 7 = 2.5 the same 
argument is true as we gradually move towards lower 
concentrations. The survival probability retains the sim- 
ple exponential character as we decrease c, even for the 
lowest trap concentrations used. The Rosenstock approx- 
imation, Eq. @, predicts this simple exponential decay 
only in the time range where (5'„) ~ n. As we have seen, 
though, in Fig. |2] there is a crossover in the behavior of 
(Sn) with time, which should modify this behavior. How- 
ever, this crossover takes place at early times and is not 
apparent in the linear time scale used for the survival 
probability. The Rosenstock approximation, thus, based 
on the results of Fig. |21 predicts a simple exponential 
decay for $ on scale-free networks. 

Fig. Efa) validates, thus, the assumption that the 
Rosenstock approximation is true in the case of scale-free 
networks with 7 = 2.5. The decay of the survival proba- 
bility, though, is greatly influenced by 7. In Fig.QIb), the 
results for 7 = 3.5 and large trap concentrations c clearly 
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FIG. 5: Survival probability as a function of time, for net- 
works of size A^ = 10® and trap concentration c = 0.01. From 
left to right, the network connectivity is (symbols) 7 = 2.0, 
2.5, 3.0 and 3.5. The corresponding Rosenstock approxima- 
tion is represented with thin lines. We also present the sur- 
vival probability for regular networks (thick lines) in (top to 
bottom) d = 1, 2 and 3 dimensions. 



demonstrate a deviation from a simple exponential be- 
havior and the failure of the Rosenstock approximation. 
Only in the case of low c, such as c = 10^"^, this ap- 
proximation is satisfactory and describes reasonably well 
the exponential decay of the simulation data. Thus, for 
7 = 3.5 we also employed the cumulant approximation 
of Eq. (jHJ. The higher-order moments of the 5'„ distri- 
bution were calculated numerically, via the same simula- 
tions that yielded the first moment (5„) of Fig. |21 It is 
evident that the description of the data improves signif- 
icantly. The second-order truncation (i.e. including the 
standard deviation of the Sn distribution) follows quite 
closely the simulation data for c = 0.05 and c = 0.005 
over more than three decades on the vertical axis. In the 
case of c = 0.01 we need to include higher moments in 
order to achieve the same level of accuracy, since $2 cap- 
tures only part of the behavior. The fourth-order trun- 
cation $4 seems to be quite succesful over almost four 
decades, and describes a significant non-exponential part 
of the curve quite well. 

In Fig.|31we present the survival probability (c — 0.01) 
in different scale-free networks and in regular lattices. We 
can see that as 7 increases the survival probability be- 
comes higher. Since the number of connections between 
the nodes decreases with 7 and we have seen that the av- 
erage value of the number of sites visited also decreases, 
the walkers will spend more time in smaller network re- 
gions. This has a dramatic influence on $ and as we can 
see in the figure the difference in the survival probability 
between networks of 7 = 2 and 7 = 3.5 can be two orders 
of magnitude, even only after a few hundred steps. The 
shape of the curves is also different, since the exponential 
character of the lowest 7 values is no longer retained for 



7 > 3. This change in the decay, along with the much 
slower relaxation is a manifestation of the network struc- 
ture, which for 7 > 3 corresponds to a loosely connected 
network where the number of nodes with extremely high 
connectivity has diminished. 

Inspection of Fig. [3 and similar simulations for differ- 
ent concentrations on networks with 7 = 3 suggest that 
in the range 2 < 7 < 3 the Rosenstock approximation 
provides a reliable description in the time regime studied 
in this work. On the contrary, when 7 > 3 this approxi- 
mation is not valid and one needs to resort to the use of 
higher moments in the cumulant expression. 

Concerning the comparison with regular lattices, it is 
obvious that trapping in the most connected networks 
(7 =2-3) behaves in a similar manner as in 3-dimensional 
lattices (simple exponential decay), and for 7 < 2.5 de- 
cays in a similar rate, too. The case of a two-dimensional 
lattice represents the borderline dimension for recurrent 
random walks in lattices, and the relaxation of $ is not 
exponential, while for d = I the survival probability is 
considerably higher, since the walkers are confined be- 
tween two trapping sites and perform a random walk in 
this region. Similarly to the d < 2 cases, the survival 
probability relaxation in networks with 7 > 3 is not ex- 
ponential and, in general, cannot be described by the 
Rosenstock approximation. 

Scale-free networks have been considered heuristically 
to behave as infinite-dimensional lattices. This assump- 
tion (d —^ 00), however, implies that both the Rosen- 
stock approximation (Eq. 0J) and the Donsker-Varadhan 
result (Eq. [TJ would yield a single exponential decay 
<i>(n) ~ exp(— n) with the number of steps n. As we 
have seen, though, this result can be verified in the pre- 
sented time scale by the simulations for networks in the 
range 2 < 7 < 3, but not for 7 > 3. The reason is that 
in d — !■ CX3 the probability for a walker to revisit a site is 
vanishingly small, since at every step the walker has an 
infinite number of possible sites to jump to. Thus, the 
revisitation probability tends to zero and the number of 
sites visited is equal to the number of steps performed 
{{Sn) ~ n). Eq. I0J then predicts the same behavior as 
O, i.e. $(n) ~ exp(— n). For scale-free networks the 
situation is somewhat different, though. Although there 
are a few highly connected nodes in the system (hubs), 
from where a walker can be directed to previously un- 
sampled areas of the network, the largest percentage of 
the nodes has a very small number of links, e.g. fc = 1 



or fc = 2. A walker that reaches such a node will return 
at the next step to its former position. The character 
of a scale- free network as a substrate for random walks, 
thus, cannot be described as purely infinite-dimensional. 
The dimensionality can be considered as a local property 
which is modified according to the area of the network 
where the walker lies in. Depending on the value of 7, 
the area sampled depends on how connected a system 
can be and how easy it is for a random walker to visit 
new nodes. For sparse networks, for example, the revis- 
itation probability increases (together with the network 
diameter) and leads to larger deviations of the above law. 



V. SUMMARY 

In this work we presented numerical results on (i?^), 
(5„), and trapping in scale- free networks, which are well- 
studied processes in many other systems. 

Mean squared displacement was found to range from 
superlinear diffusion to sublinear diffusion as we varied 7, 
while the network coverage increases almost linearly with 
time for all 7 values examined. The Rosenstock approxi- 
mation is adequate for predicting the survival probability 
in the range 2 < 7 < 3, but for higher 7 it cannot account 
for the non-exponential character of the survival proba- 
bility decay with time. In this case, we found that the 
cumulant expansion can fit quite accurately the observed 
behavior. 

The mean-field character (exhibited by the validity of 
the Rosenstock approximation) for 7 < 3 can be also 
expected in the case of these networks, since the heuristic 
arguments supporting the Donsker-Varadhan expression 
(Eq. [Tj) do not apply here. Although a walk can still 
be compact, large trap- free regions do not exist on such 
a network. The main reason is the small average path 
length between any two nodes of the structure. For any 
trap distribution, there are not any network areas where 
a walker can spend a lot of time without meeting a trap, 
since the connectivity of the network allows it to easily 
escape to a different neighborhood, where traps may exist 
in a larger local concentration. When 7 > 3, though, the 
importance of the hubs in a network diminishes, and the 
behavior resembles more that of regular low-dimension 
lattices, with prominent non-exponential decays even at 
early times. 
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